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PREFACE 

This  report  was  prepared  by  Dr.  Donald  E.  Nevel,  Research  Physical  Scientist, 

Applied  Research  Branch,  Experimental  Engineering  Division,  U.S.  Army  Cold  Regions 
Research  and  Engineering  Laboratory.  The  work  was  funded  by  the  United  Stales 
Coast  Guard. 

Technical  review  of  the  report  was  performed  by  Stephen  L.  DenHartog  and  F.  Donald 
Haynes  of  CRREL. 


ICEBREAKER  SIMULATION 


by 

Donald  E.  Ncvel 


In  designing  an  icebreaker,  we  attempt  to  choose  the  size,  shape,  weight  and  thrust  that  will  ef- 
ficiently break  ice.  In  order  to  make  this  selection  rationally,  we  must  first  look  at  the  ways  in 
which  an  ice  sheet  can  break.  Any  combination  of  these  failures  can  occur  during  actual  icebreaking. 

When  a vertical  load  is  applied  to  the  edge  of  an  ice  sheet,  the  ice  bends.  With  sufficient  load  a 
crack  develops  perpendicular  to  the  edge  of  the  ice  sheet.  Then  a circumferential  crack  forms.  Sec- 
ondary circumferential  cracks  may  also  develop  closer  to  the  load  than  the  first  one.  In  general, 
breaking  the  ice  by  bending  is  very  efficient;  the  broken  pieces  are  relatively  large. 

When  a horizontal  load  is  applied  to  the  edge  of  an  ice  sheet,  the  ice  enters  a state  of  plane  stress. 
With  sufficient  load  the  ice  cracks  perpendicular  to  the  edge  of  the  ice  sheet.  If  the  two  pieces  of 
ice  are  not  constrained,  they  are  wedged  apart.  This  type  of  icebreaking  has  been  observed  in  thin 
ice. 


If  horizontal  and  vertical  loads  are  applied  simultaneously,  the  horizontal  load  helps  bend  the  ice 
sheet.  On  the  other  hand,  we  usually  assume  that  the  vertical  load  does  not  affect  the  horizontal 
plane  stresses.  The  secondary  circumferential  cracks  in  bending  may  be  significantly  influenced  by 
these  horizontal  stresses.  These  cracks  are  important  since  they  produce  smaller,  more  reasonably 
sized  pieces  of  ice. 

Another  way  in  which  a horizontal  load  can  break  the  ice  is  by  crushing.  The  process  of  crushing 
is  not  really  understood  and  is  probably  a combination  of  tensile  and  shear  failures.  Crushing  occurs 
when  the  ship  first  hits  the  ice.  Crushing  does  no  useful  work,  and  hence  should  be  kept  to  a mini- 
mum. However,  due  to  the  shape  of  the  ship’s  bow,  some  crushing  always  occurs. 

The  effect  of  a vertical  dynamic  load  is  to  cause  a local  failure.  An  example  of  this  is  the  vertical 
oscillating  icebreaker. 

Ideally,  we  would  like  to  formulate  the  problem  mathematically  and  solve  it.  The  nearest  problem 
that  has  been  solved  is  a static  vertical  load  on  a seini-infimte  plate  on  an  elastic  foundation,  which 
was  solved  first  by  Westergaard5  and  later  by  Shapiro4  and  Nevel  J 

The  problem  of  a static  horizontal  and  vertical  load  on  a semi-infinite  plate  on  an  elastic  founda- 
tion has  not  been  solved.  For  this  case,  with  a concentrated  horizontal  force  PH,  the  differential 
equation  in  Cartesian  coordinates  is 
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where  w = deflection 

x.v  = ('artesian  coordinates 
r.O  = polar  cootdinatcs 
p = mass  density  of  water 
g = gravitational  constant 
q - applied  vertical  stiess 
84  = t:  h3l\2pg  (1  uJ) 

E = Young's  modulus  of  ice 
v = Poisson’s  ratio  of  ice 
h = thickness  of  ice. 

Striving  either  of  these  differential  equations  is  still  insufficient  since  the  boundary  conditions  during 
icebreaking  may  be  considerably  different  than  a semi-infinite  plate.  In  addition,  dynamic  effects 
have  not  been  considered  in  the  above  formulation. 

Another  closely  related  problem  is  the  infinite  wedge  on  an  elastic  foundation.  This  problem  was 
solved  by  Nevel,J  assuming  beam  theory.  Milano1  has  applied  the  results  to  the  icebreaker  problem. 
Here  again  the  boundary  conditions  of  the  wedge  are  considerably  different  than  those  occurring  dur- 
ing icebreaking. 

Rather  than  considering  the  force  required  to  break  the  ice.  White6  has  considered  the  force  that 
a ship  applies  to  the  ice.  In  this  analysis  the  ship  hits  the  ice,  crushes  it,  slides  up  on  it  and  comes  to 
rest.  The  ice  sheet  is  assumed  not  to  deflect.  The  equations  of  motion  for  the  ship  are  deiived  and 
numerically  evaluated  by  means  of  a computer.  The  results  are  discussed  in  terms  of  the  sustained 
vertical  force  and  the  extracting  thrust.  This  method  has  the  advantage  of  showing  the  influence  of 
the  ship’s  parameters.  However,  the  limitation  is  negieaing  the  deflection  and  the  breaking  of  the 
ice. 

In  contemporary  icebreakers,  the  thrust  of  the  propellers  and  the  momentum  of  the  ship  cause 
the  ship  to  ride  up  on  the  ice.  After  the  ship  has  come  to  a halt,  the  vertical  force  on  the  ice  is  due 
to  raising  of  the  center  of  gravity  of  the  ship  and  the  vertical  component  of  the  thrust.  This  vertical 
thrust  is  due  to  the  trim  of  the  ship.  The  horizontal  force  on  the  ice  is  just  the  horizontal  component 
of  the  thrust.  Although  this  horizontal  force  helps  bend  the  ice,  a vertical  force  is  more  efficient. 
Hence,  in  order  to  have  a more  efficient  icebreaker  we  intuitively  say  that  the  thrust  and  momentum 
should  be  increased,  the  stem  angle  lowered,  the  coefficient  of  friction  decreased,  and  the  crushing 
reduced.  These  intuitive  concepts  agree  with  White’s  results. 

After  the  ship  has  come  to  rest  on  the  ice,  the  thrust  is  essentially  transmitted  as  a horizontal 
force  to  the  ice.  This  could  be  transmitted  as  a vertical  force  if  the  power  were  switched  from  the 
horizontal  thrusting  propeller  to  a vertical  one.  If  the  vertical  thrusting  propeller  were  at  the  stern, 
it  is  estimated  that  one  half  of  the  thrust  could  be  exerted  vertically  on  the  ice. 

In  the  past,  bow  propellers  have  been  effective  in  breaking  the  ice.  As  far  as  I can  determine, 
they  are  effective  because  of  Bernoulli’s  principle.  That  is  to  say,  with  an  increase  in  the  water  vel- 
ocity there  is  a decrease  in  pressure.  This  causes  an  additional  downward  force  on  the  ice.  In  polar 
ice,  bow  propellers  are  not  used  because  the  ice  does  not  break  before  it  strikes  the  propellers.  A 
possible  solution  to  this  problem  is  to  place  the  propellers  inside  a tunnel  going  through  the  ship. 

This  would  pull  the  water  into  the  tunnel  near  the  bow  and  discharge  it  farther  back  along  the  side 
of  the  ship.  The  use  ol  jets  rather  than  propellers  may  be  another  consideration. 

The  icebreaking  process  is  a complex  one  which  has  not  been  adequately  analyzed.  One  method 
of  solution  is  to  simulate  the  process  by  means  of  a model.  If  all  the  important  parameters  can  be 
scaled,  then  the  size,  shape,  weight  and  thrust  of  the  model  ship  can  be  vatied  to  increase  the  effec- 
tiveness of  icebreaking.  In  particular,  one  should  determine  the  most  effective  shape  of  the  bow.  A 
large  variety  of  bow  shapes,  such  as  straight,  convex  and  concave  bows  (l-'ig.  I ) should  be  investigated. 
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Figure  1.  Bow  shapes. 


T = propeller  thrust  g - gravitational  vector 

P - force  from  the  ice  p = water  pressure 

M - mass  of  the  ship  F = water  friction  force  per  unit  area 

Figure  2.  Forces  on  an  icebreaker. 

Since  the  primary  purpose  of  an  icebreaker  is  to  break  ice,  it  seems  reasonable  to  design  the  ship 
to  do  just  this  efficiently.  Other  factors  such  as  stability,  stern  swamping,  and  extraction  thrust 
should  be  considered  after  the  efficient  bow  shape  has  been  found.  One  can  design  for  the  ramming 
condition  in  which  the  ship  comes  to  a halt,  or  for  the  continuous  condition  in  which  the  velocity 
lost  during  deceleration  just  equals  the  velocity  gained  during  acceleration. 

For  the  model  of  the  slop  one  requires  that  all  the  linear  dimensions  of  the  prototype  (full  scale 
ship)  be  scaled  by  the  same  factor.  This  factor  X is  defined  as 

/.  = X/.m  (3) 

where  /.  is  a prototype  linear  distance  and  /,m  is  a model  linear  distance.  In  this  paper  a parameter 
with  a subscript  m means  the  parameter  is  referred  to  the  model.  A parameter  with  no  subscript 
means  the  parameter  is  referred  to  the  prototype. 

In  order  to  simulate  the  slop’s  motion  we  will  consider  the  equation  describing  this.  Figure  2 
shows  the  forces  on  the  slop.  For  the  prototype  the  equation  for  the  conservation  of  linear 
momentum  is 


Mg  * T * P + f pds  + JFJs  = M ^ 

ill 


(4) 


where  V - velocity 
l - time 

its  - ship’s  differential  surface  area. 


The  linear  momentum  equation  for  the  model  ship  is 


fpm  dsm  + fFm  dsm 


M_ 


dVj, 

dt 
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For  the  water,  the  equation  for  the  conservation  of  momentum  is 


P ^ = VP+PV2v  + Pf! 


(6) 


where  n is  the  viscosity  of  water  and  D/Dt  is  the  substantial  derivative.  The  momentum  equation 
for  the  watei  of  the  model  is 
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(7) 


If  the  model  is  to  simulate  the  prototype,  the  equation  of  the  model  must  equal  the  correspond- 
ing equation  of  the  prototype  except  for  a constant  multiplication  factor.  We  find  this  factor  to  be 
unity  between  eq  6 and  7 by  comparing  the  last  terms.  Comparison  of  the  second  term  says  p/L  - 
Pm/Lm  UsinB  L=XJ-m  ,his  8'ves 


P = XPm 

Comparison  of  the  first  term  says  V/t  - Fm/rm . Using/,  = X/,m  it  yields 

- - X* 

Hence  the  velocity  is  scaled  by 

v = x'/j  F . 


(8) 


(9) 


(10) 


Comparison  of  the  third  term  shows  that  it  cannot  be  scaled.  This  term  represents  the  friction  of  the 
water  and  is  relatively  unimportant;  hence  we  neglect  it. 

Now  we  compare  the  fourth  term  and  eq  4 and  5.  Using  eq  3 and  8 we  find  that  if  eq  5 is  mul- 
tiplied by  X1  each  term  is  equal  to  the  corresponding  term  of  eq  4.  Hence  we  get 


M = X3  A/m 

m 

(11) 

f = x3  f 

m 

02) 

P = X3  F 

m 

(13) 

iu.6 
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II 

(14) 

t = v^/m 

(15) 

Since  the  friction  of  the  water  has  not  been  scaled  correctly,  we  cannot  guarantee  that  eq  14  will 
be  satisfied. 
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The  force  P from  the  ice  can  be  resolved  into  two  components  and  the  components  must  obey 
the  same  scaling  as  P of  cq  13.  We  resolve  P into  a force  /'N  normal  to  the  ship's  surface  and  a force 
P tangential  to  the  ship’s  surface.  When  the  ship  is  sliding  onto  the  ice  we  have 

= fPN  06) 

where  / is  the  coefficient  of  friction  between  the  ship  and  the  ice.  For  the  model  we  have 


f m^N 


l 'sing  /’  = X5  Pm  we  obtain 


(17) 


/ = rm 


(18) 


We  note  that  if  the  conservation  of  linear  momentum  of  the  ship  is  modeled,  the  conservation  of 
angular  momentum  will  also  be  modeled.  This  follows  since  the  angular  momentum  equals  length 
times  linear  momentum.  Of  course  the  model  is  assumed  to  have  the  same  weight  distribution  as  the 
prototype 

We  now  turn  our  attention  from  the  ship  to  the  ice.  The  force  /’from  the  ship  on  the  ice  may  be 
resolved  into  a vertical  force  Py  and  a horizontal  force  PH.  First  we  shall  consider  the  static  case  of 
the  vertical  force  only.  The  maximum  force  that  can  be  exerted  on  the  ice  depends  upon  the  bend- 
ing stress  o that  will  cause  the  ice  to  break.  The  general  solution  to  the  bending  of  an  ice  sheet  under 
a vertical  load  for  a wide  variety  of  boundary  conditions  can  be  expressed  as 


/>  \ t t t t tl 


(!</) 


where  a and  b are  the  lengths  of  the  load  distribution.  For  the  model  ice  sheet  we  have 

'Y  (20) 


U_/l_  / y m am  \n 

’ii  ’sm  'em  ’em 

' m m m m m 


m m _ 


Since  the  function  G is  the  same  function  in  eq  19  and  20,  their  values  will  be  equal  it  their  argu- 
ments are  the  same.  Thus  we  obtain 


dh 


7 

m m 


Py  (Py)m 
Utilizing  eq  1 3 we  obtain 

h = \\oJo)'/7  hm  . 

In  order  for  eq  22  to  be  valid  we  need 


= xmIVm . y/i  = yjt„ 


(21) 


(22) 

(23) 

(24) 
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*/«  = ^ m.  b/t  = bj 8m 

(25) 

b/i  = ■ 

(26) 

liquation  23  is  reasonably  satisfied  since  Poisson’s  ratio  does  not  vary  greatly  Cor  most  materials, 
and  previous  solutions  show  that  the  function  G does  not  depend  greatly  upon  u. 

Let  us  assume  that  a horizontal  distance  x in  the  ice  sheet  is  scaled  by 


* = (27» 

It  seems  desirable  to  have  0 = A,  since  the  ship  is  in  contact  with  the  ice.  However  at  this  time  we 
will  not  make  them  equal.  They  can  be  made  equal  at  a later  time  if  this  proves  to  be  convenient. 
Now  eq  24  and  25  will  be  satisfied  il 

8 = 0 8m  . (28)’ 


Using  cq  22  this  becomes 


9/ 

0 = X/8| 

C-)/4 

' 1 -v2  / 

' 0 ' 

(29) 


Using  eq  28,  we  see  that  eq  26  cannot  be  scaled.  From  previous  solutions,  it  is  known  that  b/i 
has  an  influence  on  the  value  of  G only  in  the  vicinity  of  concentrated  loads  on  thick  ice  sheets.  Al- 
though the  icebreaker  may  approximate  this  for  initial  cracking,  the  ultimate  failure  docs  not  occur 
in  the  vicinity  of  the  load,  and  hence  should  not  depend  significantly  upon  h/ 8.  Therefore  eq  26 
need  not  be  scaled. 

We  note  that  from  previous  solutions,’  the  deflection  w of  an  ice  sheet  can  be  expressed  as 


Pv  \ 8 8 8 8/ 

Since  the  arguments  of  the  function  //  have  already  been  scaled  we  obtain 

..  V2  _ »m 
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(30) 


(31) 


Using  eq  13  and  28  this  becomes 

In  order  to  consider  the  effects  of  dynamic  forces  and  the  forces  in  the  plane  of  the  plate,  we 
consider  the  complete  differential  equation,  which  is 


8«V«  + 2^  Hi? 

PH  PK  dt 1 


q : 3Jw  + 32>v  + Ny  d2w 

PK  PK  <>X2  PK  tody  PK  dv5" 


(33) 
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where  7 = density  of  ice 

Ay  = force  per  unit  length  in  the  plane  of  the  plate  in  the  x direction 
Nv  - same  but  in  they  direction 
\'xv  - shear  force  per  unit  length  in  the  plane  of  the  plate. 

There  is  a similar  equation  with  subscripts  m describing  the  modeled  ice  sheet. 

Comparison  of  the  first  terms  shows  that  if  eq  33  is  multiplied  by  02/X1  it  is  equal  to  the  modeled 
equation.  The  second  terms  yield 

P = (X3/0JV»m.  (34) 

Comparing  the  third  terms,  which  represent  the  vertical  acceleration  force  of  the  ice  sheet,  and  using 
eq  9,  22  and  32,  we  obtain 


7 = (o/om)/j  X 


v'/2 


(35) 


The  fourth  term  is  satisfied  and  the  remaining  terms,  which  represent  the  forces  in  the  plane  of  the 
plate,  yield 

P = 03Pm.  (36) 


Therefore  the  forces  in  the  plane  of  the  plate  can  only  be  scaled  correctly  when  U ~ X. 

If  0 t X,  then  the  horizontal  distances  in  the  tee  and  the  forces  in  the  plane  of  the  plate  will  not 
be  scaled  correctly.  The  effect  of  the  incorrect  horizontal  length  of  contact  area  of  the  ship  will  be 
insignificant  on  ultimate  failure.  But  the  effect  of  incorrect  horizontal  scaling  on  the  size  of  the 
broken  ice  pieces  will  be  significant.  If  the  size  of  the  broken  ice  pieces  is  not  important,  such  as  in 
the  ramming  condition,  it  is  not  necessary  for  0 = X.  However,  under  the  continuous  condition, 
where  the  size  of  the  broken  ice  pieces  determines  the  ship’s  acceleration  distance,  it  is  necessary  that 
0 = X.  In  order  that  0 = X,  we  must  choose  the  model  material  for  ice  such  that  the  following  equa- 
tion gives  a reasonable  value  of  X: 


(JWiyi 


(37) 


We  note  that  if  the  model  material  is  the  same  as  sea  ice,  then  0 > X Tor  0 * X,  the  model  will 
simulate  a prototype  that  will  break  the  ice  into  pieces  proportional  to  0xm  while  a real  icebreaker 
would  break  pieces  proportional  to  Xx m . This  means  the  model  will  have  a greater  acceleration  dis- 
tance than  the  real  icebreaker.  Furthermore  the  model  will  simulate  a prototype  that  has  honzontd 
forces  of  03Pm  while  the  real  forces  arc  X 1 l'm . Both  of  these  observations  say  that  the  model  ice- 
breaker will  do  better  than  the  corresponding  real  icebreaker  By  similar  reasoning  the  reveise  argu- 
ment could  be  made  when  0 < X. 

Saunders-Roe  used  a wax  for  the  model  ice  sheet.  The  properties  of  tins  wax  were  u ,x  m>  kgy 
cin3  and  /:  = 3.17  x 1 0 ’ kg/cm2.  The  U.S.  ( oast  (iuard  used  freshwater  ice  lor  tests  at  the  \avy 
Flectronics  laboratory  in  March  1967.  The  properties  of  this  ice  were  o 1 1 27  kg  cm2  and  /•'/ 

(I  v‘ ) = 100. S x lO’kg/cm2.  For  freshwater  ice  v - ‘/3  which  gives  / = H9.3  x I0J  kg/cm2.  Low 
salinity  ice  was  also  tried  at  N.T.I..  This  ice  showed  a greater  variation  in  thickness  than  the  freshwater 
ice.  The  properties  of  this  low  salinity  ice  were  a = I 1 .72  kg/cm2  and  /•'/(  I u2 ) 9|  .0  x I O'  kg/ 

cm2  Using  v = % we  gel  / K0.9  x 10*  kg/ cm 2 The  properties  for  sea  ice  aie  approximately  n - 

5 kg/cm2  and  /■  = 20  x l()J  kg/cm2 . 
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In  order  to  determine  which  material  is  belter  we  compute  8 as  determined  by  eq  29  for  the 
three  model  materials.  Since  Poisson’s  ratio  is  not  accurately  known  for  these  materials,  we  assume 
the  same  Poisson’s  ratio  for  all.  For  X = 36  we  get  for  wax,  freshwater  ice  and  low  salinity  ice  a 0 of 
1 10,  52.5  and  54.7  respectively.  For  X = 44  we  get  for  wax,  freshwater  ice  and  low  salinity  ice  a 0 
of  I 38.  65.8  and  68  6 respectively.  Hence  it  appears  that  freshwater  ice  will  give  the  best  simulation 
of  the  process. 

There  is  a possibility  that  various  additives  to  the  water  may  produce  an  ice  which  will  give  a 0 
value  closer  to  X However  this  remains  for  future  research. 

In  summary  the  relationships  between  a prototype  and  model  icebreaker  are  determined  by  the 
following  equations 


/.  = 

m 

(38) 

M = X’  M 

m 

(39) 

t = x3  f 

(40) 

P = X3  P 

(41) 

V = x'/j  V 

m 

(42) 

i = x'*  ,m 

(43) 

f’fm 

(44) 

h = X/2(om/o)/2  hm 

(45) 

x = 0xm 

(46) 

E 

c* 

ii 

(47) 

w = (X3/02  )wm 

(48) 

where 


In  order  to  simulate  the  vertical  acceleration  of  the  ice  sheet,  we  need 

7 = (°/Om)/j  7m  (50) 


In  order  to  simulate  the  forces  in  the  plane  of  the  ice  sheet  and  the  size  of  the  broken  ice  pieces,  we 
need  0 = X,  which  means 
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